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ABSTRACT 

In this investigation we quantify the metallicities of low mass galaxies by constructing the most 
comprehensive census to date. We use galaxies from the SDSS and DEEP2 survey and estimate 
metallicities from their optical emission lines. We also use two smaller samples from the literature 
which have metallicities determined by the direct method using the temperature sensitive [OIIIJA4363 
line. We examine the scatter in the local mass-metallicity (MZ) relation determined from ^20,000 
star-forming galaxies in the SDSS and show that it is larger at lower stellar masses, consistent with the 
theoretical scatter in the MZ relation determined from hydrodynamical simulations. We determine a 
lower limit for the scatter in metallicities of galaxies down to stellar masses of ^ l{f Mq that is only 
slightly smaller than the expected scatter inferred from the SDSS MZ relation and significantly larger 
than what is previously established in the literature. The average metallicity of star-forming galaxies 
increases with stellar mass. By examining the scatter in the SDSS MZ relation, we show that this is 
mostly due to the lowest metallicity galaxies. The population of low mass, metal-rich galaxies have 
properties which are consistent with previously identified galaxies that may be transitional objects 
between gas-rich dwarf irregulars and gas-poor dwarf spheroidals and ellipticals. 



1. INTRODUCTION 

The metallicity of a galaxy is a key physical parame- 
ter f or understanding it s history and evolution. Starting 
with iPagel et "aD (|1979f l , the practice of estimating gas- 
phase abundances from optical emission lines has allowed 
for significant advances in the study of galaxy evolution. 
Metals are produced in stars and returned to the inter- 
stellar medium (ISM) through stellar mass loss. In this 
way, the metallicity provides a record of the star for- 
mation history. However, the gas-phase abundance is 
modulated by gas flows where feedback processes can 
drive enriched gas out of galaxies and accretion of both 
primordial and enriched material can influence the gas- 
phase abundance. Observations of the relation between 
metallicity and other physical properties of galaxies and 
the evolution of these relations can place important con- 
straints on the physical processes governing galaxy evo- 
lution. 

The mass-metallicity (MZ) and luminosity-met a .llicity 
(LZ) relation was flrst observed bv lLeaueux et al.l ()1979( ) 
using nearby dwarf galaxies. Us ii ig ^5 0,000 galax- 
ies from the SDSS, iTremonti efaH (|2004f ) have estab- 
lished the relation for the local population of star- 
forming galaxies down to stellar masses of log(M*) ~ 
10*'^ Mq. The MZ relation has also been observed at 
higher redshifts (Erb ct al. 2006; Cowie & Barger 200 
Lamareille et al..i2009i : .Mannucci et al..i200& : iZahid et al 



20111 among others). From the scatter in the MZ and 
LZ relation, many authors conclude that the fundamen- 
tal correlation is between stellar mass and metallicity, 
such that lower mass galaxies have lower metallicity (e.g. 
ITremonti eraIll2004l : IZahid et al.ll2011| ). ITremonti etall 
(|2004l ) attribute this effect to stellar wind driven mass 
loss enabled by the shallow potentia l wells in low mass 
galaxies. However, iDalcanton et al.l (|2004D argue that 



observations of the local MZ relation are also consistent 
with metal-poor infall of gas. 

Several authors have suggested that gas flows alone 
cannot explain the MZ relation. Usi ng smoo thed particle 
hydrodynamical simulations, ijjrooks et al.l (J200_7.) have 
argued that in addition to mass loss driven by UV heat- 
ing and supernova feedback, a mass dependent star for- 
mation efficiency can help to explain the origin and evolu- 
tion of the MZ relation. In this scenario, due to lower star 
formation efficiencies, low mass galaxies have yet to con- 
vert much of their gas reservoir in to stars and therefore 
are less metal-rich. Alternatively, iKoppen et al.l ()2007f ) 
suggest that the MZ relation may result from variations 
of the initial mass function. They argue that in low mass 
galaxies the initial mass function does not extend to high 
stellar masses due to the low levels of star formation, thus 
resulting in fewer massive stars forming and less chemical 
enrichment. 

The metallicities of low mass galaxies (M* < 10^'^ 
Mq) may provide important clues to the origin of 
the MZ relation which may help to distinguish be- 
tween the various proposed scenarios. Many au- 
thors have inv estigated the metallic it ies of nearby 
dwarf galaxies (iSearle fc Sargent! Il972t iSkillman et al.l 
19891: iKennicutt k SkiUmanI 120011: iLee et al.l 
Hidalgo-Gamez et al.l 120031 : Ivan Zee fc HavnesI 
Zhao et al. 201^ a mong others) but few have focused 
on the MZ relation (iTainura et a l.' 2001; Lee et al.ll2006l: 
iVaduvescu et al.ll2007l ). iLee et al. (2006) extend the MZ 
relation by 2.5 orders of magnitude in stellar mass. Their 
main conclusion is that the scatter in the MZ relation is 
similar over 5 orders of magnitude. 

Recent observations of low mass, metal-rich galax- 
ies p re sent a challenge t o this conclusion (Gu et al] 
[20Cia lDellenbusch"etaII [2007L iPeeples et al.l 120081 : 
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Fig. 1. — The distribution of redshift (top left), U— B color (top 
right), absolute B-band magnitude (bottom left) and stellar mass 
(bottom right) for our selected sample of ~20,000 SDSS galaxies. 



iPetropoulou et all 1201 I f). In particular, iPeeples et al.l 
(2008) focus on a population of metal-rich dwarf galaxies 
from the SDSS which they consider outliers to the MZ 
relation. They identify these galaxies as objects nearing 
the end of their star formation activity and transition- 
ing from dwarf spirals and irregulars to dwarf spheroidals 
and elliptical s , cons istent with the scenario discussed by 
iGrebel et al.l |2003[ ). However, citing possi ble system- 
atic ov erestimates in the derived metallicities. lBerg et al.l 
(j2011h refute this result. Due to the small number of 
observations, it remains unclear whether there is a pop- 
ulation of low stellar mass, metal-rich galaxies. A popu- 
lation of galaxies, such as these, would present an impor- 
tant new test for chemical evolution models and cosmo- 
logical simulations while providing stronger constraints 
for the physical mechanisms governing the MZ relation 
and galaxy evolution. 

Here we examine the metallicities in galaxies spanning 
steUar masses from 10'^ < M^/Mg < 10^°. In Section 2 
we present the data and selection criteria used and in Sec- 
tion 3 we discuss our determination of stellar mass from 
photometry and metallicity from optical spectra. We ex- 
amine the scatter in the MZ relation determined from 
^ 20, 000 SDSS galaxies in Section 4 and in Section 5 we 
focus on the metallicities of low mass galaxies. In Section 
6 we investigate systematic uncertainties in our mass and 
metallicity determinations. We present the main results 
of this study, the scatter in the MZ relation at low stel- 
lar masses in Section 7. The physical properties of these 
galaxies are examined in Section 8 and in Section 9 we 
provide a discussion. We give a summary of our results 
in Section 10. When necessary we adopt the standard 
cosmology (iJo,^m,^A) = (70 km s^^ Mpc^^, 0.3, 0.7). 

2. DATA SAMPLE 
2.1. The SDSS Sample 



The SDSS DR7 consists of -900,000 gal axies span- 
ning a redshift range of < z < 0.7 ( Aba zauan et al.l 
l2009f) . The survey has a Petrosian magnitude limit of 
rp = 17.8 and the spectra cover the nominal range of 
3900 - 9100A. For eac h object five bands of photome- 
try, ugriz, is available (jStoughton et al.ll2002[ ). We use 
the publicly available emission line fluxes measured by 
the MPA-JHU grouiQ. We determine the local MZ rela- 
tion from — 2 0,000 galaxies using the sample selection of 
iZahid et al.l ()2011[ ). Hereafter, we refer to this relation 
as the SDSS MZ relation. 

We first distinguis h star-forming galax ies from AGN 
using the scheme of IKewlev et al.l (|2006( ). In order to 
avoid aperture effects, we require a g-band fiber aperture 
covering fraction > 30% and impose a lower redshift limit 
of 0.04 (Kewley et al. 2005). The n iedian covering frac - 
tion for the selected sample is 38%. IKewlev et al.l (|2006D 
find that the SDSS sample is incomplete at higher red- 
shifts and in order to minimize evolutionary effects we 
impose an upper limit redshift cutoff of z = 0.1. 

We require that the S/N H/3 > 3, aii23 < 2 and equiv- 
alent width of H/3 be greater than 4A. Here, cr_R23 is the 
error in the R23 parameter which is the sum of the flux 
in the [OIII]A4959, 5007 -1- [OII]A3727 divided by the H/3 
flux. We flnd that the measured relation is not strongly 
dependent on the values adopted for these cuts and we 
select in this manner to be consistent with previous stud- 
ies. In the selected sample, the median S/N of Ha and 
H/3 are 37 and 21, respectively. 

Figure [T] shows the distribution of redshift, U— B color, 
absolute B-band magnitude and stellar mass for our se- 
lected sample. This particular selection gives us a pure 
star-forming galaxy sample in a narrow redshift range. 
It covers 2 orders of magnitude in stellar mass and has 
a color distribution consistent with blue star-forming 
galaxies. 

2.2. Dwarf Galaxies Sample 

Many measurements of metallicities of dwarf galaxies 
can be found in the literature. We use the sam ples of 
Lee et al. (200 3, see references t herein'). [Gi l de Pa z et al.l 
(j2003) and .Peeples et al.l (|2008| ) . iLee et ali (,2006.) extend 
the SDSS MZ relation by 2.5 decades in stellar mass. 
Their sample is taken from the literature and con sists of 
27 dw arf galaxies. Using the models of Bell & d e Jond 
()2001| ). they de rive mass es from the B — K color and 
scale them to a iChabrieii (I2003.) IMF . In th is study we 
use the 21 galaxies in the iLee et al.l (|2006( ) sample for 
which there are metallicities determined by the direct 
method and we ad o pt th eir stellar mass determination. 

iGil de Paz et al.l (J2003") compile a sample of blue com- 
pact dwarf (BCD) galaxies. These galaxies are selected 
on basis of color, morphology and luminosity such that: 
M-B^A^fl < 1, A*s < 22 mag arcsec"^ and Mk > —21 mag. 
Here, fiB and fiR are the peak surface brightness in B 
and i?- band, respectively and Mk is the K-band magni- 
tude. Zhao et al. (2010) investigate the LZ relation for 
blue compact dwarfs using B and R-band photometry 
from this compilation with addition al K.^-band photom- 
etry obtained from several sources (iJarrett et al.l (20001 : 
INoeske et al.1l20Q3L [20051: IVaduvescu et al.ll2007D . We de- 
termine, for the first time, the stellar masses from the 

^ http://www.mpa-garching.mpg.de/SDSS/DR7/ 
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Fig. 2. — The distribution of redshift (top left), U— B color (top 
right), absolute B-band magnitude (bottom left) and stellar mass 
(bottom right ) for our s a mpl e of dwarf galaxies. For the 21 dwarf 
galaxies from I Lee et al] l)200fi ^ we do not have redshift or color 
data. 

photomet ry as described in S ection 3.2. Of the 80 galax- 
ies in the iZhao et aD ()2010[ ) sample, we select 66 with 
detections of [OIIIJA4363, allowing for metallic ities to be 
determ ined using direct met hod. Simila r to | Lee et al.l 
(|2006[ ), the metaUicities of the lZhao et~aLl ()2010D sample 
are compiled from the literature and may suffer from 
aperture effects. However, by comparing abundances 
measured fr om several d i fferen t sources in a subsample of 
the galaxies. I Zhao et all (|2010[ ) conclude that the effect is 
small. Both these samples with metallicities determined 
from the direct method reveal a population of low mass, 
metal-poor gala x ies (se e Section 5.1). 

iPeeples et al.l (|2008| ) examine the outliers from the 
mass-metallicity relation. Their sample consists of 41 
metal-rich dwarf ga laxies from the SDSS. We use the 
iPecplcs ct al.' ('2008') sample but adopt the flux values 
of Berg ct al. (2011) for the four galaxies they reexam- 
ined. Given the strict selection criteria of IPeeples et al.l 
(J2008), the "very low mass" sample of galaxies (which 
have log(MH.) < 9) only consists of 17 galaxies. In order 
to investigate possible systematic effects in the determi- 
nation of meta Uicities of low stella r mass galaxies, we 
supplement the IPeeples et al.l (|2008[ ) sample of 17 galax- 
ies with 56 low mass, metal-rich galaxies taken from the 
SDSS sample described in section 2.1. We refer to these 
56 galaxies as the supplemental SDSS sample. These 
galaxies have stellar masses <1O*'*M0 and have high 
metallicities determined by three different diagnostics, 
making them outliers from the SDSS MZ relation. We 
use several indepen dent method s of metallicity dctermi- 
nation (i.e. lYin et al. 200 7; Kobulnickv fc Kewle v 2004: 
iKewley fc Dopitall2002L see Section 3.1) in order to mit- 
igate systematic effects that may result in overestimates 
of the metallicity with any one method. 

Figure [2] shows the distribution of redshift, U— B color, 
absolute B-band magnitude and stellar mass for the 
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Fig. 3. — The distribution of the parent (grey histogram) and 
selected (black histogram) sample of galaxies from the DEEP2 sur- 
vey. The redshift distribution (top left) of the selected sample is 
roughly consistent with the parent sample. Our selection picks out 
the bluer (top right), less luminous (bottom left) and less massive 
(bottom right) galaxies as compared to the parent sample. 

dwarf galaxies sample. 

2.3. The DEEP2 Sample 

We probe significantly fainter objects than is possible 
with th e SDSS by analyzi ng galaxies from the DEEP 2 
survey (| Davis et al.l l2003| ). This allows us to probe 
deeper to find low stellar mass objects otherwise not ob- 
served in shallow broadband surveys. The majority of 
galaxies in our parent sample come from observations 
of the Extended Groth Strip. The DEEP2 team used 
the DEIMOS multi-object spectrograph on the 10m Keck 
telescope to target galaxies in 4 fields covering 3.5 square- 
degrees down to a limiting magnitude of Rab — 24.1. 
They primarily target galaxies in the redshift range of 
0.7 < z < 1.4. In order to maximize efficiency of finding 
high redshift galaxies, a color preselection approximately 
given by: B-R < 2.35x (R-I) - 0.45, R-I > 1.15, or B-R < 
0.5 was applied for galaxies in fields 2-4. Galaxies in the 
Extended Groth Strip region (field 1) were targeted with- 
out a photometric redshift color pre-selection ( Coi l et al.l 
I2OOI . 

The moderately high-resolution (R~ 5000) spectra 
have a nominal spectral coverage of 6500 — 9100A. BRI- 
band photometry for these galaxies was obtained by the 
DEEP2 team using the CFH12K camera on the 3.6 m 
Canada France Hawaii Telescope ()Coil et al.l 12004) and 
for 54% of our selected sample of galaxies ifg-band pho- 
tome try is available from the Hale telescope at Mt. Palo- 
mar (|Bundv et al.ll2006[ ). In this study, we use the third 
release of the dataQ. 

In order to estimate metallicities, we use the 
[NIIJA6584 and Ha emission lines. These lines are ob- 
served in the DEEP2 data in the redshift range of < 

2 http://deep.berkeley.edu/DR3/ 



z < 0.38. We measure the uncalibrated flux of tfiese 
emission lines by taking the area of a least-square gaus- 
sian fit. In many cases, the [Nil] line can be very weak 
(S/N < 3), so we simultaneously fit two gaussians with a 
single gaussian line-width to both emission lines. We use 
the same method to measure the [SIIJA6717, 6731 lines. 
We also measure the equival ent widths u s ing th e method 
outlined in the appendix of iZahid et al.l ()2011[ ). The er- 
ror in the fitted line flux and equivalent width comes from 
propagating the measurement uncertainties in the spec- 
trum. The median signal-to-noise (S/N) in our sample 
for Ha is 22. 

From the parent sample of 2,065 galaxies within the 
requisite redshift range, we select 777 galaxies for this 
analysis. In our selected sample, 695 are found in the 
Extended Groth Strip (field 1). By examining the pho- 
tometric imaging, we find that 51 of the sources are star- 
forming knots in larger galaxies; these sources are re- 
moved from the sample. We select galaxies that have a 
reduced-x^ of the fit less than 2 and a S/N of Ha greater 
than 3. We also require that the 68% confidence inter- 
val of our stellar mass estimate be smaller than 0.3 dex. 
These selections remove 1172 galaxies from the parent 
sample. Initially, we apply no S/N cut on the [NIIJA6584 
line in order to avoid biasing against low metallicity ob- 
jects. We remove 5 galaxies from the sample that have 
log([NII]/H a) greater than ma ximum star-formi ng ratio 
of -0.3 (e.g. lKewlevet al.ll2006[ ). Following Weiner^eLaLl 
([200?) who suggest that a substantial fraction of red 
emission line galaxies in the DEEP2 survey may be AGN 
contaminated, we remove 60 red galaxies using the blue 
and red galax y color division for DEEP2 galaxies from 
IWillmer et all (l20Q6h . 

Figure [3] shows the properties of the parent (grey his- 
togram) and selected sample (black histogram). The red- 
shift range of the parent and selected sample are com- 
parable. Similar to our selection of SDSS galaxies, our 
selection of the DEEP2 data picks out the late-type blue 
star-forming galaxies which are generally less bright and 
less massive than early-type red galaxies. The U— B, Mb 
and mass distributions in Figure [3] reflects our selection. 

3. METHODS 
3.1. Mass Determination 

We estimate galaxy stellar masses by comparing pho- 
tometry with stellar population synthesis models in order 
to determine the mass-to- light ratio which we use to scale 
the o bserved luminosity (jBell et al.ll2003t iFontana et al.l 
|2004( ). M a,Knitudes are synthe s ized fr om the stellar tem- 
plates of iBruzual fc CharlotI (|2003| ) and an IMF de- 
scribed by Chabricr (2003). The 27 models span seven 
exponentially decreasing star formation models (SFR 
oc e"*/^) with T = 0.1, 0.3, 1, 2, 3, 5, 10, 15 and 30 Gyrs 
and two metallicities (0.4 and IZ©). The stellar popu- 
lation ages rang e from to 13 Gyrs and we apply the 
extinction law of iCalzetti et al.l (|2000[ ) allowing E(B-V) 
to range from to 0.6. The median statistical error for 
the derived stellar masses, determined from propagating 
the uncertainty in the photometry, is 0.15 dex. Though 
systematic effects can be significant (jDrorv et al.l 12004 : 
iConrov eF al. 2009 ), in all the samples except for the 21 
galaxies from iLee et al.l (|2006D , we have consistently mea- 
sured the stellar masses giving us a robust relative mea- 



sure. The U and B-band abs olute magnitudes f or galax- 
ies in the SDSS, DEEP2 and I Zhao et all (|2010[ ) samples 
are synthesized using this same procedure. 

E mission lin e cont ributions are accounted for by taking 
the iKennicufrn (J1998) relation between the SFR and the 
UV luminosity which is synthesized as part of the mass 
determination procedure. In our sample, making no cor- 
rection for the emission line contributions does not al- 
ter our mass determinations significantly. We adopt the 
median of the mass distribution and take the 68% confi- 
dence interval as a measure of the error. In[Zahid ct aD 
(|2011h we compare this method with the one used by 
the MPA/JHU group to determine stellar masses of the 
SDSS galaxies. We find that the results differ by a con- 
stant offset of ~0.2 dex and that the dispersion between 
the two methods is 0.14 dex. 

3.2. Metallicity Determination 

The galaxies used in this study are not homogeneously 
observed. We therefore require several methods for deter- 
mination of metallicities which we discuss here. Further- 
more, different systematic uncertainties are associated 
with the various methods of metallicity determination, 
so when possible we use several methods for the same 
galaxy in order to assess and mitigate these uncertain- 
ties (see Section 6.2). 

Metallic i ties for t he dwarf g alaxie s sample from 
iLee et al.l (|2006D and iZhao et al.l (|2010D are determined 
using the direct method (Section 3.2.1). For the SDSS 
and DEEP2 galaxies we determine metallicit ies using the 
empirically calibrated method of I Yin et al.l ([2007, Sec- 
tion 3.2.2). To address systematic bias and uncertainties 
we apply both empirically (Section 3.2.2) and theoreti- 
cally (Section 3.2.3) calibrated methods in determining 
the scatter in the metallicities of local star- forming galax- 
ies from SDSS. We also derive empirical corrections for 
metallicities determined from the [NIIJA6584 using the 
empirically calibrated N/0 ratio (Section 3.2.4 and Sec- 
tion 6.2.2). 

3.2.1. Electron Temperature Method 

The "direct" or electron temperature (Tg) method re- 
lies on measurements of the temperature sensitive auroral 
lines (e.g. [OIII]A4363, [N II]A5755). These hues provide 
tight constraints on the electron temperature and allow 
for straightforward determination of the gas-phase abun- 
dance. The metallicitie s are calcu lated using the iterative 
scheme presented by Iz otov et al.l (200 6') . Due to temper- 
ature fluctuations and gradients, the 0+ and 0^+ ions 
reside in different zones of the HII region. Te([OIII]) is 
determined for the O^"*" ion from the [OIIIJA4363 line and 
a linear conversion is applied to determine the tempera- 
ture in the region of 0+, Te([OII]). Contributions from 
higher ionization states are considered negligible and the 
gas-phase abundance is taken as the sum of the abun- 
dance in the two zones. 

This method suffers from several drawbacks making its 
use impractical in many cases. The [OIIIJA4363 line is 
on the order of 100 times weaker than the strong oxygen 
emission lines (e.g. [OII]A3727, [OIII]A5007) observed 
in optical spectra thus necessitating high S/N measure- 
ments. In most cases, this method can only be applied 
to nearby galaxies where high S/N spectra are available. 



Due to efficient cooling of HII regions by metals, the line 
can only be observed in low metallicity HII regions (< 0.5 
Z0). Finally, in the presence of temperature gradients 
and fluctuations within HII regions, several authors have 
argued that the Tg method may lead to underesti mates of 
the metall icity (e.g. Stasihska 2002. 2005: Bresolin et al.l 
|2P06; Pch a-Guerrero et al.l 1201 1[ ). although recent work 
bv iBrcsolin et al.l ()2009() and has shown good agreement 
between direct abundances of extragalactic HII regions 
and stellar metallicit ies for blue supergiants (see also 
IKudritzkiet~an[20Tl . 

3.2.2. Empirically Calibrated Method 

Semi-empirically calibrated methods rely on calibrat- 
ing strong-line ratios to metallicities determined from 
the combination of the Te method along with detailed 
photoionization modeling. The [NIIJA6584 to Ha ratio 
is shown to be strongly correlated to gas-phase oxygen 
abundance. iPettini fc Pagell ()2004D have calibrated this 
ratio to nearby HII regions. They have parameterized 
the linear relation as 



12 -I- log(0/H) = 8.90 + 0.57 x N2 



(1) 



where N2 = log([NII]A6584/Ha). This relation is vahd 
in the range of —2.5 < N2 < —0.3 and has la intrinsic 
scatter of 0.18 dex. 

The linea r relation between N2 and metallicity param- 
eterized by IPettini fc Pagell ()2004l ) is poorly constrained 
at me taUicities below 12 + log(0/H) - 7.7. lYin et al.l 
(|2007f l place stronger constraints at lower metallicities on 
the relation between metallicity and N2 by using a sam- 
ple ~ 700 galaxies from the SDSS along with values found 
in the literature using only metallicities determined from 
the Tg method, thus providing a purely empirical cali- 
bration. They parameterize the linear relation as 



12 + log(0/H) ^ 9.263 + 0.836 x N2. 



(2) 



Their relation is valid in the range of —2.5 < N2 < —0.5 
and has a la intrinsic scatter of 0.16 dex. In this study, 
we are investigating low mass galaxies in the local uni- 
verse and therefore find it desirable to have a metallicity 
diagnostic appropriately constrained at low metallicities. 
When determining metallicitie s from iV2, we adopt the 
parameterization of I Yin et al.l ()2007D . In the appendix 
we co mpare this parameterization with IPettini fc Pagell 
(|2004f l at low stellar masses. 

3.2.3. Theoretically Calibrated Method 

A third class of metallicity diagnostics relies solely 
on calibration of strong-line ratios using photoionization 
models. These methods are not susceptible to observa- 
tional limitations imposed by empirical calibrations re- 
lying on the faint [OIIIJA4363 line. Therefore, the model 
metallicities are well constrained and the parameteriza- 
tion is well defined o yer a broad range. 

The methods o f iKobulnickv fc KewlevI (|2004D and 
iKewley fc Dopital ([2002') calibrate strong-line ratios to 
metallicities determine d using photoionization models. 
IKobulnickv &: KewlevI (2004) rely on the R23 diagnos- 
tic, which is the sum of the flux of [OIIJA3727, 3729 and 
[QUI] A4959, 5007 normalized to the H/3 flux. Though 
IKobulnickv fc KewlevI ()2004t ) apply explicit corrections 
for the ionization parameter variations, there still may 



be some systematic e ffects that are unco rrected. There- 
fore, we also use the iKewlev fc Dopital (2002) method 
which relies on the iV202 ratio which is the [NIIJA6584 
to [OIIJA3727, 3729 fiux ratio. This diagnostic is shown 
to be independent of the ionization parameter. Fi- 
nally, for the SPSS s ample of galaxies we also use the 
iTremonti et al.l ()2004[) method which relies on fitting the- 
oretical models of integrated galaxy spectra to the most 
prominent optical lines. For a review of these methods, 
we refer the reader to the appendix of IKewlev fc EllisonI 
12008). We provide conversions between two commonly 
used i?23 diagnostics and the direct method in the ap- 
pendix of this paper. 

3.2.4. Nitrogen to Oxygen Ratio 

In order to investigate the effect of enhanced nitro- 
gen enrichment with respect to oxygen, we measure the 
nitrogen abundance from emission line fluxes. N/0 is 
empirically calibrated using star-forming galaxies in the 
SPSS with di r ect m easurements of ionic abundances by 
lAmorin et al.l ()2010D . The calibration is given by 

log(N/0) = -0.86 -I- 1.94 X ^^25*2 + 0.55 x N2S2'^, (3) 

where 

,.^^0 , f [NII1A6584 \ 

N2S2 = log ( roT^J^i^ onn. ] ■ (4) 



V[SII]A6717,673iy ' 

[SII] is strongly correlated to [OH] and due to less uncer- 
tainty in the reddening correction, the relation between 
N/0 and [Nil] /[SII] shows less dispersion (^ 0.1 dex) 
than the relation between N/0 and [Nil] /[OH]. 

3.3. Ionization Parameter 

The ionization parameter characterizes the level of ion- 
ization of the gas in an HII region. The ionization pa- 
rameter is defined as 



F, 



(5) 



where Fi is the fiux of ionizing photons per unit area and 
n is the number density of hydrogen atoms. The ioniza- 
tion parameter has units of velocity and can be thought 
of as the maximum velocity attainable by an ionization 
front driven by the local radiation field. The ioniza- 
tion parameter can best be constrained observationally 
by comparing two ionization states of single species. A 
commonly used diagnostic for measuring the ionization 
parameter is the 032 parameter defined as 



^o„ , /[OIII]A4959,5007\ 
^'2 = ^°S l [OII]A3727,3729 J 



(6) 



However, 'Kewl ev fc Popital (2002) have shown that the 
[OIII]/[OII] ratio is not only dependent on the ionization 
parameter, but is also sensitive to metallicity. 

Conversely, many of the strong line diagnostics are sen- 
sitive to the ionization parameter. One approach to dis- 
entangling the relationship between metallicity and the 
ionization parameter is to assign an initial guess of the 
metallicity in order to constrain the ionization parame- 
ter. Using this estimate of the ionization parameter, a 
second, more accurate estimate of the metallicity can be 
obtained. This process can be iterated until the values of 
both the ionization parameter and metallicity converge. 
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Fig. 4. — The mass metallicity relation for SDSS galaxies. The 
black squares are median metallicity in 50 bins of stellar mass de- 
termined from ~20,000 galaxies in the SDSS. The dashed line is a 
fit to the linear portion of the MZ relation (log(M*) < 9.4). The 
solid and dotted gray curves are the 68 and 95 percent contours of 
the distribution. 

We me asure the ionization pa r amete r and metallicity 
using the iKobulnickv fc KewlevI (|2004D diagnostic. This 
method uses the iterative scheme described above. The 
ionization parameter is determined from the 032 line 
ratio and the metallicity is determined from the i?23 line 
ratio, defined as 



i?23 = log(12mi^5007 



H/J 



[OII]A3727,3729 \ 



(7) 

For details of how metallicity and i onization parame- 
ter are determined see the appendix of iKewlev fc EllisonI 
(|2008h . 

3.4. Star Formation Rates 

The star formation rates (SFRs) for the SDSS samples 
were made available in the D R7. These SFRs are mea - 
sured using the technique of iBrinchmann et al.l (I2004D 
who e mpirically correct for aperture effects. iSalim et al.l 
(|2007f ) found that in galaxies with low level star forma- 
tion, the SFRs we re overestimated using this technique. 
ISalim et al.l (|2007f ) correct this overestimate and this has 
been applied to the S FRs provided in th e DR7. 

For the DEEP2 andlZhao et al.l (I2010D samples we use 
the conversion of ' Kennicuttl (|1998[ ) to determine SFRs 
from Ha luminosities. This is given by 

SFR (Mo yr~i) = 4.6 x lO"'^^ Lho (ergs s'^), (8) 

where Lhq is the Ha luminosity and the pre-factor 
has been s c aled d own by a facto r 1.7 t o convert from 
a Sal ^eteJ (ITOSSl) to a IChabrieri (l2003h IMF. For the 
DEEP2 data we scale th e uncalibrat e d Ha fluxes to the 
photometry and for the iZhao et alj (120101) samp l e, Ha 
luminosities are determined bv iGil de Paz et al.l ()2003l ) 
from narrow-band Ha filter images. We correct for dust 
extinction in both these samples using the E(B— V) de- 
rived from the SED fit when determining stellar mass. 

4. THE SDSS MZ RELATION AND ITS SCATTER 

We determine the SDSS MZ relation from --20,000 
galaxies with metallicities determined using the A'^2 



method given by Equation [51 In Figure JH we show 
the MZ relation and its scatter. We note that through- 
out this paper, when referring to the scatter in metal- 
licities, we explicitly mean to refer to the range of ob- 
served metallicities at a fixed stellar mass. The relation 
is determined by sorting the data into 50 equally pop- 
ulated bins of stellar mass and taking the median mass 
and metallicity in each bin. The relation is plotted by 
the black squares and the scatter is shown by the solid 
and dotted grey curves which are the 68 and 95% con- 
tours of the distribution, respectively. For galaxies at a 
fixed stellar masses < lO^'^M© the scatter is symmet- 
ric. At higher stellar masses, the distribution is skewed 
towards lower metallicity galaxies (with a skewness near 
— 1). A large part of this skewness can be attributed to 
the saturation of stroiig-line methods at high metallici- 
ties (|Kewlev fc Dopital l2002f). Additional skewness may 
also be the result of a physical upper limit to the metal- 
licity attainable by a galaxy, though this is not well es- 
tablished and separating out these two effects is beyond 
the scope of this work. 

The dashed line in Figure JH is a linear fit to the data. 
The turnover observed in Figure 0] is due to a combina- 
tion of a real turnover in metallicities as a function of 
stellar mass and saturation of the N2 parameter at high 
metallicities. We only fit to the linear portion of the 
SDSS MZ relation below the turnover (12 + log(0/H) < 
8.8). The fit is performed by minimizing the square of 
the residuals and is parameterized as 

12 -I- log(0/H) = (8.585±0.003)-l-(0.47±0.01)XAf. (9) 

Here, Xm = log{M^) — 9 and Af* is the stellar mass in 
solar mass units. We fix the zero-point at a stellar mass 
of 1O^M0 in order to reduce the covariance between the 
slope and intercept. The errors in the fit are assessed 
by bootstrapping the sample. We note that we consider 
this fitted MZ relation as the fiducial relation for later 
c omparisons with the data. 

iLee et all ()2006[ ) examine the MZ relation for 27 low 
mass galaxies and conclude that the scatter is simi- 
lar over five decades of stellar mass. However, using 
^ 20, 000 galaxies from SDSS (FigureH]), it becomes clear 
that the scatter increases substantially for lower stellar 
mass galaxies down to stellar mass of '^10^'^ M©. We 
test whether the larger scatter at lower stellar masses 
is real or an artifact of the method used in determining 
metallicities. We examine the scatter as a function of 
stellar mass using several diagnostic methods of metal- 
licity determination. We quantify the scatter in the MZ 
relation as the magnitude of the interval containing 95% 
of the data for a given stellar mass bin. 

Figure [5] shows the sc atter as a function of s tellar mass 
for several diagnostics. iTremonti et al.l ()2004| ) determine 
the MZ relation from ~50,000 galaxies in the SDSS. 
The solid black line is the scatter in their determina- 
tion of the MZ relation. We have subtracted 0.2 dex 
from their mass determination for consistency w ith our 
measurements of stellar mass (Zahi d et al.l 120111 ). The 
green dot-dashed, red triple dot-dashed and blue dot- 
ted curves are t he sc atter for the MZ relation, with 
the I Zahid et alj (120111) selection, determined using the 
diagno stics of [K obulnicky & Kcwlc)^ ()2004[ ). lYin et ahl 
(|2007f) and [Kewlev fc Dopitai ([20020, respectively. The 
blue long-dashed curve is the theoretical scatter de- 
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Fig. 5. — The width of the metalhcity interval containing 95% 
of the data in each mass bin as a function of stellar mass for 
different metallic ity diagnostics. The solid black line is the 
data taken from 'Tremonti et al.l l|2004l ). The dot-dashed green 
curve, triple dot-dashed red curve and the dotted blue curve 
are the Zahid et al. (2011) selected data with motallicities deter- 
mined using the Kobulnicky & Kewley (2004), Yin et al. (200'3) 
and IKewlev fc Dopita (.2002 ) diagnostics, respectively. The blue 
long-dashed curve is the theoretical scatter from hydrodynamical 
simula tions using a momentum conserving wind model (|Dave et al.i 
12011^ 1. 

termined from cosm oloKJcal hydrodynamical simulations 
(jDave et al.ll2011aL more details in Section 9). 

Strong line methods, in particular diagnostics using 
A^2, are known to saturate at high metallicities. This 
is partly responsible for the low scatter observed at 
high stellar masses for the MZ relations determined 
from the N2 diagnostic (red curve). However, the 
trend of decreasing scatter in the MZ relation is con- 
tinuous and is observed across all stellar masses, in 
particular at lower stellar masses where the strong- 
line diagnostics are not saturated. The mean rela- 
tive errors in the N2, N202 and R23 line ratios used 



in the lYin et al.l (120071) . IKewlev fc Dopital |2002) and 
iKobulnickv fc KewlevI (|2004[ ) diagnostics are 0.02, 0.05 
and 0.04 dex, respectively. The observational uncertain- 
ties of the N202 and i?23 line ratios increase with stellar 
mass due to the diminishing line strength of the oxygen 
lines at higher metallicity. Even at high stellar masses, 
the observational uncertainties are substantially smaller 
than the observed scatter and we note that for SDSS 
galaxies, systematic uncertainties associated with strong- 
line methods dominate over the observational uncertain- 
ties (jKewlev fc Ellison! [200 81 . In general, the magnitude 
of the observed scatter in Figure [S] is significantly larger 
than either the observational or systematic uncertainties 
and we attribute this to intrinsic scatter in metallicities 
of galaxies at a fixed stellar mass. 

In Figure m we see that the reason the scatter is larger 
at lower stellar mass is due to the fact that the metallicity 
of the most enriched galaxies in any given stellar mass 
bin decreases with a much more shallow slope than the 
least enriched galaxies. When comparing the high mass 
end to the low mass end of the MZ relation in Figure IH 
we see that the metallicities of the most enriched galaxies 
only decrease by ~0.2 dex over 2 decades of stellar mass, 
whereas the metallicities of the least enriched galaxies 
decrease by ~1 decade. 
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Fig. 6. — The MZ relation for 21 dwarf irre gular galaxies(red 
triangles) from the sample of I Lee et aLl (I2006|) and 66 c ompa ct 
blue galaxies (green stars) from the sample of lZhao et al.l l(20l3). 

5. METALLICITIES OF LOW MASS GALAXIES 

In this section we extend our study of the scatter in 
the MZ relation to low mass galaxies. Because different 
methods have been applied in determining the metallic- 
ities of metal-poor and metal-rich galaxies, we present 
them separately in Sections 5.1 and 5.2, respectively. 

5.1. Metal-Poor Dwarf Galaxies 

We plot the mass a nd metallic it y of th e sam ple of dwarf 
galaxi es taken from iLee et ahl ()2006[ ) and I Zhao et all 
()2010[) in Figure m The metallicities for this sample have 
been determined using the direct method (see Section 
3.2.1). The fiducial MZ relation is plotted as the dashed 
line. In Figure [6] we linearly extrapolate the 68 (solid 
line) and 95% (dotted line) contours to lower masses by 
fitting the contours in the linear portion of the fiducial 
relation (8.5 < log(M*) <9.4). We extrapolate the ob- 
served scatter to investigate whether at the low mass 
end the scatter is consistent with the expectation from 
the fiducial relation. 

The slope of the MZ relation determined from this sam- 
ple alone is 0.27 ± 0.02. The fit and error have been 
determined by a bootstrapping method. The slope dif- 
fers significantly (9cr) from the slope of the fiducial rela- 
tion (0.47 ± 0.01, see Equation JH]) determined from the 
N2 method as shown in Figure JH We note that the 
N2 method used in determining the fiducial relation is 
calibrated to be on the same absolute scale as the di- 
rect method. The metallicities of the 87 dwarf galaxies 
plotted in Figure [6] deviate strongly from the distribu- 
tion inferred from the SDSS MZ relation such that no 
low mass, metal-rich galaxies are observed. The smaller 
sc atter observed in low stellar mass galaxies was noted 
bv lLeeet all (|200l . 

The small scatter observed in these data could be a re- 
sult of selection bias attributable to the direct method of 
metallicity determination. The [OIIIJA4363 line becomes 
too weak to observe at '^0.5 Z© (12 + log(0/H)^8.5) and 
above. This is because the [OIIIJA4363 line strength is 
anti-correlated to metallicity, such that lower metallicity 
objects have stronger emission. These observational ef- 
fects could lead to an artificial suppression of the scatter. 
In particular, the fact that all the galaxies in the sample 



have 12 + log(0/H) < 8.5 and that the highest metal- 
hcity galaxy observed at a given stellar mass is a strong 
function of stellar mass likely results from observational 
biases associated with the direct method. 

5.2. Metal-Rich Galaxies from SDSS and DEEP2 

In Figure [7] we plot the masses and metallicities 
of the DEE P 2 galaxies along with the sample from 
IPeeples eFal.. (,2008^ and our supplemental SDSS sample. 
The metallicities for these galaxies ar e determined from 
the A'^2 calibration of lYin et al.l (|2Q07[ ) given in Equation 
H For 31% (240/770) of the DEEP2 sample of galax- 
ies [N1IJA6584 is not detected at a ia level. For these 
galaxies we adopt a icr upper limit value for [NIIJA6584 
when determining the metallicity. The dashed line is the 
extension of the fit to the linear portion of the SDSS 
MZ relation as seen in Figure [S) The 68 (solid line) and 
95% (dotted line) contours of the fiducial relation have 
been extrapolated for stellar masses below log(M*) < 
8.5. The median observational uncertainty in the metal- 
lici ties of the detected DEEP 2 galaxies is 0.05 and for 
the IPeeples et al.l ()2008| ) and SDSS supplemental sample 
it is 0.008 dex. 

Opposite to the behavior of the [0111JA4363 line, the 
[N11JA6584 line becomes stronger at higher metallicities. 
The fixed magnitude limit of the DEEP2 survey means 
that at lower stellar masses only the strongest [NIIJA6584 
emitters will be be detected. This effect can be seen in 
Figure[7]where the lowest metallicity object observed at a 
given mass increases in metallicity with decreasing stellar 
mass. This is op posite to t he observa tional bias present 
in the samples of JLee et all (p006) and lZhao et al.l (|20100 
which have direct method metallicities. 

These data for the first time reveal a population of 
low mass, metal-rich galaxies down to stellar masses of 
~10^ Mq. By combining two methods in determining 
metallicity which have the opposite observational bias, 
we find a significantly higher scatter in the metallicities 
of low mass galaxies than previously reported. Because 
of incompleteness and sample bias, we cannot reliably 
determine the MZ relation or quantify its scatter at low 
stellar masses. In Section 7 we derive a lower limit for the 
observed scatter in the metallicities of low stellar mass 
galaxies. 

6. SYSTEMATIC UNCERTAINTIES IN MASS AND 
METALLICITY 

6.1. Stellar Mass Estimates 

Our apparently low mass, metal-rich objects could be 
higher mass galaxies with systematically underestimated 
stellar masses. This could possibly result from heavy ob- 
scuration due to dust. For our sample of 94 metal-rich 
galaxies from SDSS, we determine the extinction from 
the Balmer decrement. The mean E(B— V) for our sam- 
ple is 0.29 ± 0.12. We do not observe H/3 in our DEEP2 
sample and therefore determine the extinction from the 
SED fit. The mean E(B-V) for our DEEP2 sample is 
0.26±0.16. Our metal-rich sample of galaxies have values 
of extinction consi stent with star-form ing galaxies in the 
local universe fe.g. lJansen et al.ll200lD and are not found 
to suffer from heavy extinction. Moreover, our method 
for determining stellar masses corrects for extinction and 
we remove galaxies from our sample that have stellar 



mass estimates that have 68% confidence intervals > 0.3 
dex. 

The velocity dispersion of a galaxy has been shown to 
correlate with the stellar mass with some scatter (e.g. 
iKassin et al.ll2007n . Therefore, low velocity dispersion 
can be taken as an indication of low stellar mass. The 
line-of-si ght velocity dispersi on for the DEEP2 sample is 
given bv lWeiner et al.l (|2006l ) as 



^obs 



(10) 



where (Jdisp is the velocity dispersion given in km s^^, 
Xobs is the observed wavelength and Uobs and Gmst are 
the observed line width and instrumental resolution, re- 
spectively. For the DEEP2 sample, the velocity disper- 
sion is primarily determined from the Ha emission line. 
For the DEEP2 data the instrumental resolution is 0.56 A 
(jWeiner et al.ll200"6l ). The velocity dispersion becomes 
unreliable when (Jdisp > cJinst/^obs due to the fact that 
small errors in the line width translate to large errors in 
the velocity dispersion (Wciner et al. 2006). 

Figure [5] shows that the stellar mass is correlated to the 
velocity dispersion, as expected. The median velocity 
dispersion of our DEEP2 sample is 26 km s~^. Below 
^^25 km s~^, shown by the vertical dashed line, galaxies 
have unreliable determination of the velocity dispersion 
due to their small line widths but are still consistent with 
low stellar masses. We are not able to perform the same 
comparison for our sample of galaxies from SDSS because 
of the lower instrumental resolution of that survey. 

6.2. Metallicity Estimates 

Recently, iBerg et ahl (|201lD have reexa mined four of 
the metal-rich galaxies from the sample of IPeeples et al.l 
(2008i 'l. Using new spectroscopic observations along 
with detail e d com parisons with photoionization models, 
IBerg et ahl ()201lD conclude that the metallicities of the 
four galaxies have been overestimated. They attribute 
this overestimate of the oxygen abundance to high N/0 
and to low ionization. 

6.2.1. Ttie Ionization Parameter 

We gain a handle on the ionization level of the gas 
by examining the ionization parameter (see Section 3.3 
for derivation) . In Figure |9] we have plotted the ioniza- 
tion parameter as a function of ste llar ma s s. Th e red 
points plot the sample from iPeeple ^ et al.l (l2008f) with 
the four galaxies inve stigated bvlBerg et al. l (|2011l ) cir- 
cled. For half of the IPeeples et al.l ()2008l ) sample, the 
[Oil] doublet is not observed. For these data we have 
used the [SII]A6717, 6731 line as a proxy The [SII] flux 
is strongly correlated to the [Oil] flux owing to their sim- 
ilar ionization potential and primary origin. We derive 
a linear relation between the [SII] an d [Oil] dereddened 
hue flux from -20,000 galaxies in the lZahid erall ()201lD 
sample and use this to infer the [Oil] line flux when only 
[SII] is observed. The rms in the relation is 0.16 dex. 
These objects are shown by the red stars in Figure El 
The dashed black curve is the ionization parameter for 
the SDSS sample binned by stellar mass and the gray 
solid and dotted curves are the 68 and 95% contours, 
respectively. 

We are not able to ascertain the distribution of ioniza- 
tion parameter for our DEEP2 sample due to inability 
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Fig. 7. — The mass and metallicity for our sample of metal-rich galaxies. The black points are 534 galaxies from the DEEP2 survey. 
The gray arrows are an additional 248 galaxies from DEEP2 where we only measure upper limits. These upper limits arc determined by 
adopting a 3a upper limit for non-detected [NIIJA6584. The red and blue squares are galaxies from the SDSS taken from Poeples et J] 
11200 8) and our supplemental SDSS sample, respectively. The median error in the metallicities of the detected DEEP2 and SDSS sample of 
galaxies are 0.05 and 0.008 dex and are shown by the black and red error bars in the top left, respectively. The dashed line is the fit to the 
SDSS MZ relation take from Figure (6] The solid and dotted gray lines are the 68 and 95% contours of the SDSS MZ relation. The fit and 
contours have been extended down to lower masses by linearly extrapolating the contour from the linear portion of the SDSS MZ relation 
(log(M.) < 9.4). 

galaxies have low ionization, though this is not a general 
feature of low mass, metal-rich galaxies. Future observa- 
tions may shed light on this issue.^ 

Consistent with the findings of iBerg et al.l ()2011[ ). we 
find that three of the four galaxies that they investigated 
have ionization parameters that are lower than the larger 
SDSS sample. However, many of the galaxies in the sup- 
plemental SP S S sam ple (blue points) and the sample of 
iPeeples et al.l ()2008l red points) appear to have ioniza- 
tion parameters consist ent with the larger SDSS sample 
from lZahid et al.l ()2011[ ). This suggests that though low 
ionization may be an issue for some of the low mass, 
metal-rich galaxies, it is not generally the case. 

6.2.2. Empirical Correction for Nitrogen Enhancement 

For our metal-rich sample we use the [NIIJA6584 line 
to measure the oxygen abundance. Nitrogen is known 
to have a primary component, formed mostly in massive 
stars, and a secondary component, formed in low and 
intermediate mass stars. Oxygen on the other hand is a 
primary element. Oxygen abundances inferred from A'^2 
line ratio are known to depend on nitrogen to oxygen ra- 



2 


: ' 


' 


= 


1 1 


^ 




i 


10 


L 


• '• • 


J 


9 


L . ..{^ 


<.^^* '. • 


: 


8 


.'•'■ ■^'i^y-t^^a^^- 


A 








7 


h • • . • „.:,i', 


, 


i 


6 


r 




\ 


1=1 


\ , 





: 



Velocity Dispersion (km/s) 

Fig. 8. — The stellar mass plotted against the velocity dispersion 
for our sample of DEEP2 galaxies. The velocity dispersion below 25 
km/s (dashed line) has large errors due to instrumental resolution. 

to measure the same line at two ionization stages (e.g. 
[Oil] and [OIII]). It may be that our sample of DEEP2 
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Fi g. 9. — The ionizatio n parameter plotted against stellar mass. 
The lPeeples et al.l l l200l) sample for which [OII]A3727, 3729 is ob- 
served is plotted by the red squares. When [Oil] is not observed, 
we use [SII] as a proxy. These galaxies are plotted by red stars. 
The 56 galaxies comprising our supplemental sample are plotted 
by the blue squar es. The 032 value for SDSS sample taken from 
IZahid et al.l 1)20111 ') is plotted by the dashed line. The data is sorted 
into 50 bins of stellar mass and the median 032 value is taken in 
each bin. The 68 and 95% contours of the data arc plotted by the 
solid and dotted gray curves. 



tip (N/0) CStorchi-Bergmann ct al."1994l :lDenicol6 et all 
I2002t (Pere z-Montcro & Contini 200^^ The secondary 
production of nitrogen is dependent on the amount of 
oxygen aheady present in the star through the CNO 
cycle, producing a larger N/0 at higher metallicity. 
Similarly, this dependence leads to a low dispersion in 
the N/0 ratio at low metallicities (jEdmunds fc Pagell 
119781 : lAUoin et al.l 11979 1) with the dispersion increasing 
at higher metallicities |Pcrez-Montero & Contini 2009). 

We determine N/0 using Equation O For the DEEP2 
sample we do not have flux calibrated data, so we sub- 
stitute in the equivalent widths of [Nil] and [SI I ] lines 
whe n determining 7V25 '2. iKobulnickv fc Phillips! l)2003D 
and IZahid et"an ()2011| ) have demonstrated that substi- 
tuting equivalent widths for line fluxes when measuring 
line ratios does not introduce significant errors. Figure 
[TU] shows the N/0 plotted as a function of stellar m ass. 
The nitrogen abundances of the lPeeples et al.l (|2008[ ) and 
supplemental SDSS samples are elevated relative to the 
DEEP2 sample. This is consi stent wi th the generally 
higher metallicities of the Peeples et al.l |2008) and sup- 
plemental SDSS samples as compared to the DEEP2 
sample (see Figure [71 and due to metallicity depen- 
dency of secondary production, an independent measure 
of these galaxies having high metallicities. 

We derive an empirical correction for enhanced ni- 
trogen abundance by comparing the difference between 
direct method metallicities with those measured using 
the strong line method. We select all star-forming 
galaxies in the SDSS sample that have S/N > 5 de- 
tections of [OII]A3727,3729, [OIII]A4363, 4959, 5007, H/3 
and [SIIJA6717, 6731. There are 627 galaxies that meet 
this criteria. We measure th e oxygen a b undan ce using 
the strong line calibration of lYn et al.l (l2007f). th e di- 
rect method as parameterized W 



N/0 using the calibration of lAmorin et al 

Figure [TT] we plot the difference in metallicity between 

the strong line method and the direct method as a func- 
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Fig. 10. — N/O plotted as a function of stellar mass. The black 
points are 373 galaxies from the DEEP2 sample. For 59 galaxies 
in the DEEP2 sample, we observe [SII] with S/N < 3. For these 
galaxies we have adopted a Str limit for the [SII] EW. These data 
give a lower limit for N/O and are plott ed by the gra y arrows. 
The red and blue squares are the i Peeples et al.l II2008I) and our 
supplemental sample, r espectively. The f our galaxies circled are 
the ones reexamined bv IBerg et al.l 120111) . The median error for 
the nitrogen abundance of 0.1 and 0.05 dex for the DEEP2 and 
SDSS sample are shown by the black and red error bar in the 
bottom right corner, respectively. 
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Fig. 11. — The difference in metallicity between the strong line 
method and the direct method plotted against N/O. The sample 
is 627 galaxies from the SDSS DR7 with emission lines of inter- 
est detected with a S/N > 5. The red line is a linear fit to the 
data taking into account errors in both coordinates. The median 
error in the metallicity difference and N/O is 0.09 and 0.03 dex, 
respectively. 



tion of N/O. The strong line method increasingly overes- 
timates the metallicity with respect to the direct method 
as a function of N/O for this sample. 

We have determined the errors in Alog(0/H) by 
adding in quadrature the errors from the strong line 
and direct method metallicities. The median error in 
Alog(0/H) and N/O is 0.09 and 0.03 dex, respectively. 
The direct method metallicities dominate the errors in 
Alog(0/H) owing to the weakness of the [OIII]A4363 
line. We fit a linear relationship to the data using 
the MPFITEXY routine which takes into acco unt er- 
rors in both coordinates ([Williams et al.l [201 Cft. The 
MPFITEXY routine depends on the MPFIT package 
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(|Markwardtll2009f ) . The linear fit is given by 

Alog(0/H) = (0.64 ± 0.03) + (0.42 ± 0.02) x log(N/0) 

(11) 
and is plotted by the red line in Figure [H] The rms of 
the fit to the data is 0.11 dex, with 0.07 dex attributable 
to the observational uncertainties. 

6.2.3. Correction Overestimates 

In making the correction for N/0, we have assumed 
that the direct method provides the most reliable metal- 
licity estimate and therefore have corrected our strong 
line metallicities. We offer some words of caution in 
making this type of correction. The correction has been 
derived as a function of N/0. The galaxies used in 
deriving the empirical correction mostly lie in the re- 
gion between -2.0 < log(N/0) < -1.0 (see Figure [T^ . 
We have extrapolated the correction for galaxies with 
log(N/0) > —1.0, w here most of the galaxies from the 
iPeeples et al.l (|2008l ) and our supplemental SDSS sample 
lie. 

In applying this correction, we have assumed that 
the direct method provides the most reliable metallicity 
measure. However, temperature fluctuations and gra- 
dients in HII regions may lead to underestimates of the 
metallicity with the direct method fStasiiiska 2002, 2005; 
iBresolinet 311 12006'). The [OIII]A4363 line strength in- 
creases with temperature. In the presence of fluctuations 
or gradients the inferred temperature may be biased to- 
wards higher temperatures in which case it would not 
be repre sentative of the HII reg ions, especially in global 
spectra (jKobulnickv et al.lll999| l. 

7. SCATTER IN MZ RELATION AT LOW MASSES 

In Figure [T^] we apply our derived empirical correc- 
tion for enhanced nitrogen enrichment to the data. Even 
with the correction for enhanced N/0, low mass, metal- 
rich galaxies remain. The empirical correction for N/0 
bri ngs the metalliciti es of the four galaxies investigated 
by IBerg et al.l (l20TTh down to 12 + log(0/H)~8.4-8.5 
(circled galaxies in Figure [T^ . The low ionization in 
three of the four g alaxies wou l d furth er reduce the metal- 
licity estimate. IBerg et all (|2011[ ) conclude that the 
metallicities of their four galaxies probably lie in the 
range of 7.9 < 12 -I- log(OH) < 8.4. Our analysis supports 
this conclusion. However, the same analysis applied to 
the rest of the galaxies in this study also supports the 
conclusion that not all of the low mass, metal-rich galax- 
ies have overestim,ated metallicities. 

For galaxies above 10^ Mq, Figure fT2l suggests that we 
may be overcorrecting for N/O. Because of this potential 
overcorrection and the fact that the sample is incomplete, 
the scatter at low stellar masses could only be larger than 
what we observe. We consider Figure [T^] to display a 
lower limit to the intrinsic scatter in metallicities as a 
function of stellar mass. The observed scatter down to 
lO'' Mq is only slightly lower than what is inferred from 
the SDSS MZ relation. 

8. PHYSICAL PROPERTIES 

We explore some of the physical properties of the low 
mass galaxies and compare them with typical galaxies in 
the SDSS sample to understand their physical nature. 



8.1. Galaxy Colors, SFRs and Equivalent Widths 

The color-magnitude diagram is shown in Figure 1131 
The low mass galaxies have lower luminosities than the 
main sample of SDSS galaxies (bla ck contours) and tend 
to be slightly bluer in color. The IPeeples et all (|2008D 
galaxies (red squares) appear to be the reddest galaxies of 
all the samples. The metal-rich galaxies from the DEEP2 
survey and the supplemental sample from SDSS appear 
to have color s that are consiste nt with the galaxies from 
the sample of lZhao et al.l ()2010f) , which are blue compact 
dwarf galaxies selected on the basis of their color. The 
gray contours are for ^ 200, 000 galaxies from SDSS se- 
lected to be non-AGN using the BPT diagram and have 
z < 0.1. The red sequence is shown by the locus of galax- 
ies in_the_to2riglrt_of the figure. Some of the galaxies 
of IPeeples et al.l (|2008| ) have colors consistent with the 
red sequence, supporting the conclusion that these are 
transitionary objects. 

In Figure [T4l we examine the SFRs as a function of stel- 
lar mass. High mass galaxies tend to have higher SFRs. 
The SFRs of Peeples et a l. (2008) sample generally are 
lower than the other samples, consistent with their red- 
der colors. The distribution of SFRs for the DEEP2, 
I Zhao et all ()2010f ) and supplemental SDSS samples ap- 
pear to be similar. The black dashed curve is the SFRs 
for ~ 140,000 star-forming galaxies in the SDSS sorted 
into 100 stellar mass bins. These data are selected to 
have a S/N>5 in Ha and H/3, z < 0.1 and are required 
to be classified as star-forming in the BPT diagram. In 
general, the data used in this study appeare to be biased 
towards higher SFRs relative to the local galaxies from 
SDSS. For the DEEP2 data, this is partly a resuh of red- 
shift evolution in the SFRs as these data have z < 0.4. 
This plot illustrates some of the bias in the samples used 
in this study. A more complete sample is required to 
probe galaxies with lower SFRs. 

In Figure [Tsj we plot the metallicities of our galaxies 
against the equivalent width of Ha (EW Ha). Here we 
have not applied our correction for elevated N/O as this 
may be an overestimate and does not affect the inter- 
pretation. The EW Ha is a measure of the Ha line fiux 
divided by the underlying continuum. The flux of Ha 
is a proxy of star formation and the underlying contin- 
uum is dominated by the emission of the older, low and 
intermediate mass stars. This population probed by the 
underlying continuum of Ha dominates the stellar mass 
of a galaxy. The EW Ha represents a measure of the 
amount of star formation normalized to stellar mass and 
therefore changes with time. We choose to plot EW Ha 
instead of specific star formation rate because it is a di- 
rectly observable quantity and is not subject to the un- 
certainties in stellar mass and star formation rate. There 
appears to be a trend in each of the samples separately 
such that galaxies with lower EW Ha tend to have higher 
metallicities. 

8.2. Physical Nature of Low Mass Galaxies 

The blue colors observe d the low mass galaxies are 
typical of dwarf galaxies. iSchombert et al.l ([1995) find 
that their sample of dwarf spiral galaxies have colors that 
are bluer than the normal early-t ype spirals despite thei r 
low star formation rates (see also lHidalgo- Gamez"2004') . 
These colors are also consistent with those observed in 
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Fig. 12. — The inetallicity plotted against stellar mass. For the DEEP2 (black dots and gray arrows VPeepl es et al.l IJ200SI . red squares) 
and supplemental (blue squares) sample we have determined met allic ities using the 7V2 diagnostic. We apply an e mpirical correction 
to these da ta for enhanced nitrogen enrichment given in Equation 1111 We also plot the samples of ILee et aLl II2006I . red triangles) and 
IZhao et al.l l|20iQ', green stars) for which metallicities have been determined using the direct method. 
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Fig. 13. — The col or-magnitude diag ram of our sample of galax- 
ies. The DEEP2 and I Zhao et al.l 120101) s ample are plot t ed in black 
dots and green stars, respectively. The IPeeples et aLl 120081 ) and 
supplemental low mass, metal-rich sample from SDSS are plot- 
ted by the red and blue square, respectively. The black contours 
are for the main S DSS sample of star-forming galaxies taken from 
IZahid et aEI 1I20T1I') . The gray contours for ~ 200,000 galaxies in 
the SDSS. 



quiescent gas - rich d warf galaxies (|van Zee et al.l [l997l ). 
IPeeples et al.l (|2008[ ) sample of galaxies are, in general, 
the reddest galaxies. This is consistent with the lower 
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Fig. 14. — The SFR p l otted as a function of stellar mass. The 
DEEP2 and IZhao et al.l 120101) sample are plotted in black dots 
and green stars, respectively. The Peoples et al. (2008) and sup- 
plemental low mass, metal-rich sample from SDSS are plotted by 
the red and blue square, respectively. The black dashed curve is 
the median SFR in 100 bins of stellar mass for ~ 140, 000 star- 
forming galaxies in the SDSS. The 68 and 95% contours are shown 
by the solid an dotted gray curves, respectively 

SFRs observed in those galaxies. 

Blue colors are a general property observed in star- 
forming galaxies. Blue colors imply a population of 
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Fig. 15. — The metallicity as a function of the eq uivalent width of 
Ha (EW Ho). The DEEP2 and lZhao et~all 1 120101) sample are plot- 
tedin black dots and green stars, respectively. The Peoples et a\] 
l|2008l ) and the supplemental low mass, metal-rich sample from 
SDSS are plotted by the red and blue square, respectively. The 
median metallicity in 50 bins of EW Ha for the SDSS sample taken 
fromFZahid et al. (2011) is plotted by the dashed black curve. The 
68 and 95% contours are shown by the solid an dotted gray curves. 
The equivalent width is observed to strongly decrease with increas- 
ing metallicity. 

young, high mass stars which tend to dominate the lu- 
minosity output of a galaxy. The fact that blue colors 
are observed in these dwarf galaxies even though many 
exhibit low levels of star formation as inferred from the 
Ha line luminosity suggests that even low levels of lo- 
calized star formatio n in low stellar mas s systems can 
dominate the color. iPeeples et al.l (|2008[ ) olaserve that 
many of their galaxies are dominated by blue cores. In 
the supplemental SDSS sample, a few galaxies also have 
blue cores, though in the majority of galaxies star for- 
mation appears to be more widespread. 

One of the most interesting features of our sample of 
galaxies is the relation between metallicity and EW Ha 
as seen in Figure [T 5} A similar trend is observed in the 
DEEP2, SDSS and lZhao et~all (12010) sample such that 
the highest metallicity objects in each sample have the 
lowest EWs. However, though the three samples have a 
similar range in EWs, they are somewhat segregated in 
metallicity such that they cluster into what appears to 
be distinct populations. In Figure [TSl there appears to 
be a continuum between the DEEP2 and SDSS data, but 
a gap in metallicity between the Zhao et al. (2010) sam- 
ple. We hypothesize that a complete census of low mass 
galaxies cover the full range of observed metallicities and 
EWHa. 

The low and high metallicity sample of galaxies inves- 
tigated in this study have similar masses, colors, SFRs 
and EW Ha. This implies that these galaxies likely have 
similar stellar populations. Under this assumption, the 
most straight-forward interpretation of Figure [15] would 
be that the differing metallicities in the samples is due to 
the gas content of the galaxies, such that the lower metal- 
licity galaxies at a given EW Ha are more gas ric h. This 
is con sisten t with the inte r preta tion of iDellenbusch et al.l 
(|2007f) and IPeeples eFail (|200l . 

9. DISCUSSION 

Understanding and quantifying the scatter in the MZ 
relation is crucial to uncovering the origin of the MZ 



relation. Observations of low mass, metal-rich galaxies 
help shed light on this issue but also present us with a re- 
cently discovered population of galaxies requiring further 
investigation and characterization. 

The MZ relation at low stellar masses has only been 
studi e d by a few autho r s (Tamura et al. 2001c iLee et al.l 
120061: iVaduvescu et al.l 120071: iPetropoulou et al.l 120 iC 
These studies have been hampered by a small number 
of observations owing to the f aintness of low st ellar mass 
galaxies and selection biases. iLee et aP (|2006t ) conclude 
that MZ relation extends to lower metallicities with no 
increase in the scatter. We have shown that the scat- 
ter in the MZ relation increases substantially (~0.4 dex) 
over the 2 decades of stellar mass observed in the SDSS 
MZ relation (Figure [5|). 

By examining data from the DEEP2 and SDSS sur- 
veys, we are able to identify low mass, metal-rich galax- 
ies that appear to be consistent with the increasing scat- 
ter implied by the SDSS MZ relation (Figure i]). The 
observations of low mass, metal-rich galaxies challenge 
the notion that the MZ relation extends to low stel- 
lar ma sses with no increase in the scatter as found by 
ILee et al. (2006). Furthermore, given their low stellar 
masses and shallow potential wells implied by the veloc- 
ity dispersion (see Figure [S]), it is difficult to understand 
these galaxies in the context of the canonical explanation 
for the MZ relation which argues that shallow potential 
wells in lower mass galaxies prevent them from enriching 
due to loss of metals through galactic scale outflows (e.g. 
lTremontieraII[200l . 

The MZ relation is observationally well established 
in both the local (e.g. iTremonti et al.l I2004D and high 
redshift universe (e.g. lErb et al.l l2006t iMannucci et al.l 
2009). The MZ relation shows that the median (or mean) 
metallicity of star-forming galaxies increases as a func- 
tion of stellar mass. Our results suggest that in the local 
universe this is largely driven by a more rapid decline in 
the metallicity of the most metal-poor galaxies relative 
to metal-rich galaxies at lower stellar masses (see Figure 
[7] or I12p . This observed trend puts stronger constraints 
on theoretical chemical evolution models attempting to 
reproduce the MZ relation. 

In Figure [5l we have compared the scatter in the MZ re- 
lation from several diagnostics with the theoretical scat- 
ter taken from th e smoothed pa rti cle hydrodynamica l 
simulations of iDave et al.l (|2011al lbl). iDave et alJ (|2011aD 
investigate the MZ relation using several prescriptions for 
winds. Their simulations include a no wind model, two 
constant wind models (with velocities of 340 and 680 km 
s~^) and a momentum driven wind such that the wind 
speed scales with the velocity dispersion of the galaxy. 
The constant wind models have considerably higher scat- 
ter at low stellar masses than the observations and the no 
wind model generally shows less scatter than the obser- 
vations. The momentum driven wind model best charac- 
terizes the observational data and the scatter predicted 
by this model is shown in Figure [S] The overall normal- 
ization is a free parameter in the models but in addition 
to the scat ter, the shape and slo pe of the local MZ rela- 
tion (fr om iTremonti et al.ll200^ are also well fit by this 
model (|Dave et al.ll2011ar see their Figure 1). The pre- 
dicted scatter in the models deviates most strongly at 
lower stellar masses. Therefore, the scatter in the MZ 
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relation at lower masses will provide important tests for 
theoretical models of chemical evolutio n. However, cur- 
rent simulations (e.g. iDave et 3111201 lallbl ) do not resolve 
galaxies with stellar masses below 10^ M© and we await 
higher resolution simulati ons to make this co mparison. 

From their simulations. IDave et al.l ()2011aD argue that 
galaxies evolve around an equilibrium MZ relation. Large 
scale gas flows perturb galaxies from equilibrium and the 
scatter in the MZ relation is a quantitative measure of 
the time it takes to return to equilibrium which, in their 
simulations, is reestablished through infall of gas. They 
argue if the dilution time of the gas (gas mass divided by 
inflow rate) is small compared to the dynamical time of 
the galaxy (at the virial radius), then the scatter will be 
small. Within this framework, our observations suggest 
that the timescales for galaxies to equilibrate once they 
are perturbed are longer at lower stellar masses (i.e. the 
dilution time is longer than the dynamical time at lower 
stellar masses). 

In the local universe, the SFRs of star-forming galaxies 
are (anti-)correlated to the metallicity in the sense that 
the higher metallic ity objects at a given stellar mass have 
lower SFRs (Ellison et all 120081: iMannucci ct al. 2010; 
IDave et al.ll2011al ). For our low mass samples, this trend 
is only weakly present in the data (see Figure [HI). Specif- 
ically, the scatter in the MZ relation at low stellar masses 
can be slightly reduc ed if SFRs a re accounted for in the 
manner presented in IMannucci e t al. (2010). However, 
the SDSS supplemental sample and some of the DEEP2 
galaxies appear to have both high metallicities and SFRs, 
inconsistent with the observed trend at higher stellar 
masses. The full implication of this correlation is still 
not well understood, though it is probably related to the 
gas content which regulates both the gas-phase abun- 
dance and the SFR. The fact that this effect is relatively 
weak in lower mass galaxies suggests that the relative 
contribution of various physical mechanisms governing 
the metallicities of low mass galaxies may differ from 
higher mass galaxies. 

Our analysis supports the conclusion that the sample 
of low mass, metal-rich galaxies appear to be "transi- 
tional" dwarf galaxies nearing the end of their s tar for- 
mation ([GrebdiraD [IMl EeeilepEaD [Ipol) . This 
interpretation suggests that metal-rich galaxies have low 
gas fractions relative to other dwarf galaxies of similar 
stellar mass and the low gas fractions are largely respon- 
sible for the high gas-ph ase oxygen abundance observed. 
iDellenbusch et aLr(|2007f) come to similar conclusions for 
their smaller sample of low luminosity, metal-rich galax- 
ies. In this evolutionary scenario, these galaxies would 
represent the link between the gas-rich dwarf irregulars 
and the gas-deficient dwarf spheroidals and ellipticals. 
The bluer galaxies in the supplemental SDSS sample 
would be at an earlier evolutionary stage in this tran- 
sition, where s tar fo rmation is more widespread than the 
iPeeples et al.l (|2008[ ) sample but gas content is lower than 
their low metallicity counterparts. Accurate determina- 
tions of the gas content of low mass, metal-rich galaxies 
are important to test and establish this scenario. 

The environment of these low mass galaxies may play 
a crucial role in determining their chemical properties. 
Both observations and simulations have suggested that 
galaxies in dense environments tend to have metallic- 



ities over and above those expected fro m the correla- 
tion of stellar mass and the environment (jCooper et al.l 
l2008t [Dave et al]|2011al) . Moreover, the stripping away 
of gas from low mass galaxies in high density environ- 
ments, either through ram pressure stripping or inter- 
actions, has been cite d as the possible explanation for 
high metallicities (Bos elli et al. 120081 : iPetropoulou et al.l 
1201 It ). However. , Peeples et al.l (|2008D conclude that their 
galaxies are fairly isolated and the DEEP2 and supple- 
mental SDSS sample of galaxies, at least in projection, 
do not appear to be in high density environments though 
further investigation into the environments of low mass, 
metal-rich galaxies is necessary. 

In order for ram pressure stripping to occur effec- 
tively it is estimated that Pram '^ Pigmv'^ > cr^pgas/S 
(jGunn fc GottlfTOT^ . Here, Pram is the ram pressure, 
PiGM is the gas density of the intergalactic medium 
(IGM), V is the velocity of the galaxy through the IGM, 
a is the velocity dispersion and pgas is the gas density of 
the galaxy. Typical IGM densities (nn < 10~^ cm~^) are 
likely not su fficient to remove si gnificant amounts of gas. 
However, as lGrebel et al.l (|2003l ) argue, a combination of 
tidal effects along with ram pressure stripping as dwarf 
galaxies pass near giant galaxies or dense regions of the 
IGM (if the medium is inhomogeneous) is the most plau- 
sible mechanism for gas removal. In this scenario, even 
if dwarf galaxies are not currently found in high den- 
sity environments, previous encounters could remove gas 
from the ISM and any subsequent star formation would 
lead to elevated levels of enrichment due to the low gas 
fractions. 

The blue cores and low SFRs of many of the metal-rich 
galaxies suggests that star formation is not widespread 
throughout the galaxy. Localized star formation opens 
up the possibility that the high metallicities observed 
in some of these galaxies may be a result of the neb- 
ular emission dominated by individual HII regions. In 
this case, the inferred metallicity may not be ind i cative 
of the global metallicity of the galaxy. iGu et ahl (|2006[ ) 
were the first to observe a metal-rich dwarf galaxy. From 
detailed photometric and spectral analysis of IC 225, 
they find that there are two photometric cores in this 
galaxy that are spatially separated; the off-nuclear core 
being bluer than the blue nucleus. From the redshifts of 
the line emission and absorption, they conclude that the 
line emission originates from the off-nuclear core. Follow 
up kinematic study of the nuclear region of IC 225 us- 
ing integral field spect roscopy have confirmed this result 
(|Miller fc Rudiell2008l) . 

In most star-forming galaxies, the gas-phase oxygen 
abundance is determined from emission lines integrated 
over many star-forming regions and hence considered to 
reflect the global metallicity. If, in these low mass sys- 
tems, the emission lines are dominated by a single region, 
the metallicity determined from these lines may not re- 
flect the global gas-phase abundance. Spatially resolved 
spectroscopic studies of the metal-rich dwarf galaxies will 
help to establish whether the high metallicities observed 
are representative of global metallicities or are localized 
to individual HII regions. If the latter is true, determin- 
ing at what stellar masses this effect becomes pronounced 
and the systematic effects this has on metallicity studies 
of galaxies is crucial. 
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10. SUMMARY 

In this contribution, we have investigated the metahic- 
ities of low mass galaxies. Our data conies from culling 
the literature and from identifying low mass galaxies in 
the SDSS and DEEP2 survey. We summarize our main 
findings as follows: 

1. We examine the scatter in the local stellar 
mass- met allicity relation determined from ~20,000 
galaxies in the SDSS. We find that the scatter 
increases with decreasing stellar mass. We ob- 
serve this trend using several diagnostic methods 
of metallicity determination mitigating the sys- 
tematic effects and uncertainties associated with 
any particular diagnostic. The theoretical scatter 
taken from smoothed particle hydrodynamical sim- 
ulations with momentum driven winds is also ex- 
amined and found to be consistent with the ob- 
served scatter (see Figure [5]). 

2. The low stellar mass MZ relation derived from 
galaxies with metallicities determined using the di- 
rect method shows a smaller range of metallicities 
than the SDSS MZ relation. We attribute this 
lower observed scatter to selection bias in galaxies 
with metallicities determined by the direct method. 
Using the N2 strong-line diagnostic method, we 
provide observations of low mass, metal-rich galax- 
ies. We examine systematic uncertainties of en- 
hanced nitrogen enrichment and variations in the 
ionization parameter to help establish the metal- 
richness of these low mass galaxies. The results 
challenge the notion that the scatter in the low 
stellar mass-metallicity relation is small and con- 
stant with stellar mass. Due to incompleteness of 
our sample and the possibility that our empirical 
correction for elevated N/0 is overcorrecting the 
data, we give a lower limit to the scatter at low 
stellar masses (see Figure [T2|). 

3. Low mass, metal-rich objects have been identified 
in the literature as the transitional objects be- 
tween gas rich dwarf irreg ulars and gas poor dwarf 
spheroidals and ellip ticals (jDellenbusch et all2007l : 
iPeeples et al.l 120081 ). We examine the physical 
properties of these low mass galaxies and find that 
they are generally consistent with this scenario. 

Further observations are required to understand these 
recently discovered population of galaxies. In particular, 
accurate gas mass determinations will help to establish 
these low mass, metal-rich galaxies as gas poor. Spatially 
resolved spectroscopy of individual galaxies will help to 
determine if the metallicities inferred from the emission 
lines are indicative of the global metallicities. The ef- 
fects that nitrogen enhancement and low ionization pa- 
rameter have on metallicities are not well understood. A 
firm observational basis is required for deriving empirical 
corrections. Finally, detailed investigations into the en- 
vironments of these galaxies may shed light on physical 
mechanisms leading to their ostensibly high metallicities 
and their roles within the sequence of galaxy evolution. 
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APPENDIX 
Al. Comparison of N2 Diagnostics 

In Figure [16} w e compare the diagnost ics of 



iPettini fc Pagell (|2004 left panels) and lYin et aD (|2007l 
right panels). In the top panels we compare data com- 
piled from the literature tak en from the objects inves- 
tigated by iLee et all (|2006D . The compiled data are 
emission line measurements for individual HII regions. 
Only HII regions where the [OIII]A4363, [NII]A6584 and 
Ha lines are measured are included. In the bottom 
panels, the comparison between the methods is made 
for galaxies taken from the samp les of ILee et al.l ()2006f ) 
and (Zhao et~all ()2010D . For the ILee et al.l (|2006[) sam- 
ple, the integrated flux is determined by summing the 
flux from individual HII regions. The Te metallicities 
were calcula, t ed us ing the iterative scheme presented by 
llzotov et aP (|2006f ). 

In Figure [jni the solid line in each plot is the one- 
to-one agreement. The data are independent and have 
errors in both the x and y direction. We perform a linear 
bisector fit using the routine robust Jinefit. pro in IDL in 
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Fig. 16. — A comparison of the N2 metallicity diagnostic with 
the direct method determination using the [OIII]A4363 temper- 
ature sensitive line (Izotov et al. 2003)- We have determined 
the metallicities using the strong line diagnostics calibrated by 
[Pettini & Pagel (2004, left) and Yin et al. (2007, right). The top 
panels show metallicities derived for individual HII regions. The 
bottom panels show global metallicities . The red triangles are 20 
dwarf irregular galaxies from I Lee et all (|2006l ) . For these galaxies 
we sum individual spectra for a given galaxy to obtain an inte- 
grated measure. The gr een stars are 60 com pact blue galaxies with 
metallicities taken from [Zhao et all (|2010l ') and the N2 values are 
taken from their paper. The solid line in each plot is the one-to-one 
agreement and the dashed line is a linear bisector fit. 

order to assess how well the data derived from the diag- 
nostics agree with the direct metallicity determination. 
The dashed lines are a linear bisector fit in each case. 
A scatter between the relation of the direct and strong- 
line method of 0.30, 0.35, 0.23 and 0.24 dex are observed 
in the top left, top right, bottom left and bo ttom right 
panels , respectively. The greater slope in the lYin et al. 



(2007?) calibration as compared to the iPettini fc Page! 



( 2004) calibration leads to a slightly lar ger scatter (com- 
pare E quations [T] and [2])- Though the IPettini fc Pagell 
(|2004[ ) metallicit ies match the dir ect determinations rea- 
sonably well, the lYin et al.l ()2007| ) metallicities provide a 
more consistent calibration at lower metallicities. 



A2. Conversion of R23 to Te Metallicities 

We give the conversions required to make metallicities 
from two commonly used calibrations of i?23 consistent 
with the direct method at low metallicities. The left and 
right panel of Figure [IT] show the metallicities for the 
same HII r egions as F i gure 1161 de termined from the cal- 
ibratio ns of lMcGaugiil (|199lD and iKobulnickv fc KewlevI 
(|2004f ). respectively. The relatively low metallicities of 
these HII regions place them on the lower metallicity 



branch of the i?23 diagnostic. T he slope of the h near 
bisector fit is 0.67 and 0.57 f or the iMcGaughl (|199lD and 
IKobulnickv fc KewlevI (|2004[ ) calibrations respectively. 

At low metallicities R23 diagnostics present some ad- 
vantages over N2. The reduced scatter in the i?23 metal- 
licity calibration as compared to the N2 calibration is 
due to observational uncertainty. At lower metallici- 
ties, the oxygen lines are significantly stronger than the 
[NIIJA6584 line, whose line strength diminishes with de- 
creasing metallicity. This places some limitations on the 
usefulness of the A^2 diagnostic for estimating metallici- 
ties in high redshift galaxies where we expect less enrich- 
ment and therefore lower metallicities. 



9.0 



8.5 



M91 /^'' 
* •JtfSf* ** * 

' ^ ^ *x * 


KK04 y^' 

' ' ' ■ / 



7.0 



7.5 8.0 

12 + log(0/H) 



7.0 



7.5 8.0 8.5 

12 + log(0/H), T, 



9.0 



Fig. 17. — The metallicities determined from i?23 compared with 
the direct method. The l eft an d right pan el show the metallicity 
determ ined using the .McGaugK (1991!) and IKobulnickv &: KewlevI 
H2004I ) metallicity diagnostics, respectively. The solid line is the 
one-to-one correspondence and the dashed line is a linear bisector 
fit to the data. 



In order to make the metallicity found for the lower 
branch of the i?23 diagnostic consistent with the direct 
method for our sample of low metallicity HII regions, we 
provide the following conversion: 



M91' = 1.60 X M91 - 4.87 
KK04' = 1.90 X KK04 - 7.70. 



(12) 
(13) 



Here, M91 and KK04 are t he metallic i ties d eter- 
mined using the diagnostics of IMcGaughl ()199lD and 
IKobulnickv fc KewlevI ()2004[ ). respectively. The prime 
quantities are those diagnostics linearly converted to be 
consistent with the direct method metallicities deter- 
mined from the iterative scheme of llzotov et al.l (J2006D . 
This is also consistent with the A^2 line ratio metallicity 
calibration determined by lYin et aLJ ()2007D . We note 
that these conversions should only be used when the 
metallicities are known to be on the lower metallicity 
branch of i?23. 
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